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ABSTRACT 

We explore the mass-assembly and chemical enrichment histories of star forming galax- 
ies by applying a population synthesis method to a sample of 84828 galaxies from the 
Sloan Digital Sky Survey Data Release 5. Our method decomposes the entire observed 
spectrum in terms of a sum of simple stellar populations spanning a wide range of ages 
and metallicities, thus allowing the reconstruction of galaxy histories. A comparative 
study of galaxy evolution is presented, where galaxies are grouped onto bins of neb- 
ular abundances or mass. We find that galaxies whose warm interstellar medium is 
poor in heavy elements are slow in forming stars. Their stellar metallicities also rise 
slowly with time, reaching their current values (Z* ~ 1/3Zq) in the last ~ 100 Myr 
of evolution. Systems with metal rich nebulae, on the other hand, assembled most 
of their mass and completed their chemical evolution long ago, reaching ~ Zq 
already at lookback times of several Gyr. These same trends, which are ultimately 
a consequence of galaxy downsizing, appear when galaxies are grouped according to 
their stellar mass. The reconstruction of galaxy histories to this level of detail out of 
integrated spectra offers promising prospects in the field of galaxy evolution theories. 
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1 INTRODUCTION 

One of the major challenges of modern astrophysics is to 
understand the physical processes involved in galaxy forma- 
tion and evolution. Significant steps in this direction could 
be made by tracing the build up of stellar mass and metal- 
licity as a function of cosmic time. 

One way to address this issue is through cosmologi- 
cally deep surveys which map how galaxy properties change 
for samples at different redshifts (z). Among these proper- 
ties, the relation first observed by Lequeux etal. (1979) be- 
tween heavy-element nebular abundance and galaxy mass, 
or its extension, the luminosity-metallicity relation (e.g., 
Skillman etal. 1989; Zaritsky etal. 1994), are being exten- 
sively used to probe the metal enrichment along cosmic 
history. Clear signs of evolution are being revealed by 
studies of these relations at both intermediate (Savaglio 
et al. 2005; Lamareille et al. 2006; Mouhcine et al. 2006) and 
high z (Shapley et al. 2005; Maier et al. 2006; Erb et al. 2006), 
which generally find significant offsets in these relations 
when compared to their versions in the local Universe. 

Alternatively, one may study galaxy evolution through 
fossil methods, which reconstruct the star formation his- 



tory (SFH) of a galaxy from the properties of their stars. 
The most detailed studies of this sort are those based on 
nearby galaxies (including our own), where color-magnitude 
diagrams and even spectroscopy of individual stars are fea- 
sible (e.g., Smecker-Hane etal. 1996; Hernandez etal. 2000; 
Dolphin 2002; Rizzi et al. 2003; Skillman et al. 2003; Monaco 
et al. 2005; Koch etal. 2005). Comparing such observations 
with predictions of evolutionary synthesis models then al- 
lows one to infer how stars were born as a function of time, 
as well as to constrain the chemical enrichment histories of 
the galaxies (e.g., Lanfranchi, Matteucci & Cescutti 2006; 
Carigi, Colin & Peimbert 2006). By mapping how stars be- 
come richer in metals as a galaxy ages, these studies provide 
useful constraints for chemical evolution models. 

Beyond the local group and its neighborhood, fos- 
sil methods must derive SFHs from integrated galaxy 
spectra, using either selected spectral features (Kauff- 
mann et al. 2004; Brinchmann et al. 2004; Gallazzi et al. 2005; 
Nelan et al. 2005; Thomas et al. 2005; Bernardi et al. 2006) or 
the full spectrum (Panter, Heavens & Jimenez 2003; Heav- 
ens et al. 2004; Cid Fernandes et al. 2005; Mathis, Chariot & 
Brinchmann 2006; Panter etal. 2006). This approach works 
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with much less information than is available from individ- 
ual stars in galaxies, so one cannot expect similar levels of 
accuracy. On the other hand, integrated spectra of galax- 
ies are now available by the hundreds of thousands from 
modern spectroscopic surveys, so the decrease in accuracy 
is compensated by orders of magnitude increase in statistics 
and by the broader range in galaxy types spanned by such 
surveys. It is therefore of great interest to explore methods 
to retrieve information on galaxy evolution from integrated 
spectroscopic data. 

In our ongoing series of papers entitled the Semi Em- 
pirical Analysis of Sloan Digital Sky Survey (SDSS) Galax- 
ies, (Cid Fernandes et al. 2005; Mateus et al. 2006a; Stasihska 
et al. 2006; Mateus et al. 2006b; SEAGal I-IV, respectively) 
we have shown that a decomposition of galaxy spectra in 
terms of instantaneous bursts of different ages (t*) and 
metallicities produces both excellent spectral fits and 
astrophysically meaningful results. So far, however, our de- 
scription of SFHs has been extremely sketchy, based on av- 
erages over the t* and distributions, thus washing out 
valuable time-dependent information. 

This paper shows that we can actually do much better. 
To illustrate this we have culled a sample of star-forming 
(SF) galaxies from the SDSS (Section|5J, to which we apply 
our synthesis method and derive mass-assembly and chemi- 
cal enrichment histories (Section 01. The evolution of stellar 
mass and metallicity is then studied by grouping galaxies 
with similar nebular abundances or mass (SectionQJ. Details 
of this analysis and cross-checks on the results reported here 
are presented in a companion paper (Asari et al. , in prep, 
hereafter A07). 



2 DATA 

We have recently completed detailed pixel-by-pixel spectral 
fits of 573141 galaxies from the Main Galaxy Sample of the 
SDSS Data Release 5 (Adelman-McCarthy et al. 2006) with 
the synthesis code STARLIGHT described in SEAGal I and 
II. STARLIGHT decomposes an observed spectrum in terms 
of a sum of simple stellar populations covering a grid of 
25 ages from t+j = 1 Myr to 18 Gyr and 6 metallicities 
Z* ]3 - = 0.0001 to 0.05 (ie, from 0.005 to 2.5 Z e ). Each of 
these TV* = 150 instantaneous bursts is represented by a 
spectrum extracted from the evolutionary synthesis models 
of Bruzual & Chariot (2003) for a Chabrier (2003) initial 
mass function, Padova 1994 tracks, and STELIB library (Le 
Borgne et al. 2003). Illustrative fits are presented in SEAGal 
I and A07. 

Emission lines were measured from the starlight sub- 
tracted spectra. Figure 01 shows the Baldwin, Phillips & 
Terlevich (1981, BPT) diagnostic diagram for 175234 galax- 
ies for which H/3, [OIII]A5007, Ha and [NII]A6584 are all 
detected with S/N > 3. Reddening corrected line fluxes are 
used in this plot and throughout our analysis. The distri- 
bution of points in this diagram resembles a flying seagull, 
with SF galaxies occupying the left wing and galaxies host- 
ing an AGN on the right. We define SF galaxies as those 
under the divisory line proposed in SEAGal III (solid line 
in Figure^,). Further requiring a minimum continuum S/N 
of 10 at ~ 4750 A leaves a set of 84828 objects, hereafter 
the SF-sample, shown in Figure 03. Results for alternative 



sample definitions (involving different S/N thresholds, other 
AGN/SF separation criteria, and cuts on galaxy inclination) 
are presented in A07. Note that by selecting SF galaxies we 
indirectly select spirals over elliptical galaxies. 

The SF-wing is essentially a sequence in nebular metal- 
licity (SEAGal III and references therein), which we quantify 
by the oxygen abundance obtained through the strong-line 
calibration of Stasihska (2006): 

where (O/H) = 4.9 x 10~ 4 (Allende Prieto, Lambert & As- 
plund 2001). We believe this is the most reliable calibration 
available in the literature, since it is based on temperature- 
based abundances in the metal-poor case and on the ([Oil] 
+ [OIII])/H/3 calibration of Pilyugin (2000) in the metal-rich 
regime, using a large sample of giant HII regions. 

Z ne f> is both a physically motivated and mathematically 
convenient coordinate to map galaxy positions along the left 
wing in the BPT diagram. From the tip of the SF-wing in Fig 
01 to its bottom, our SF sample spans the Z ne b ~ 0.2-1.6 
Zq range. In Fig 03 this interval is chopped onto 6 loga- 
rithmic bins of width A log Z ne b = 0.13 dex, except for the 
first one which is twice as wide to include more sources. 
Galaxies inside these same bins will be grouped together 
in the analysis of star-formation and chemical enrichment 
histories presented in the following. 



3 METHOD OF ANALYSIS 

As thoroughly discussed in the literature on stellar popula- 
tion synthesis, the combination of noise in the data, astro- 
physical and mathematical degeneracies limits the amount 
of information about galactic histories that can be extracted 
from integrated spectra. Clearly, the 150 components of the 
population vector used in our fits are a greatly over-detailed 
description, so some a posteriori compression of the results 
to a lower dimension is needed. Our previous papers (SEA- 
Gal I-IV) have taken this compression approach to its ex- 
treme by reducing the whole t+ and Z t distributions encoded 
in the population vector to their first moments, 
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where /ij is the fraction of the total stellar mass (M*) which 
is in the j population. 

Both (logt*) and {Z t } proved to be mathematically ro- 
bust and astrophysically valuable summaries of the full pop- 
ulation vector. In SEAGal I, II and A07 these averages are 
shown to cleanly reveal important relations, such as those 
between (logt*) and M* (associated with the galaxy down- 
sizing phenomenon), (Z+) and M* (the mass-metallicity re- 
lation), and (Z*) and Z ne b (which reveals a connection be- 
tween stellar and nebular chemical enrichment levels) . Given 
this success, it is natural to explore less drastic compressions 
of the population vector. 
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log [NII]X6584/Ha 

Figure 1. (a) 175234 SDSS galaxies in the BPT diagram. The solid line (Stasihska et al 2006) divides Star-Forming galaxies from those 
hosting AGN. The dashed line shows the dividing line used by Kauffmann etal. (2003). (b) The SF sample, chopped into six bins of 
nebular abundance. The number of galaxies in each bin is given in the right, along with the corresponding mean Z ne b values (in solar 
units). Galaxies close to bin borders are not plotted for clarity. 



One way to compress the population vector but retain 
the age-information is to marginalize over the Z^j distribu- 
tion and compute the total mass converted into stars as a 
function of time: 

where ^ denotes the fraction of the total mass converted 
onto stars (M+) that is associated with the j base age. 
The superscript "c" stands for "converted onto stars" , and is 
introduced to distinguish fij and M* from fij and M*, which 
refer to the mass still locked inside stars, ie, discounting the 
mass returned to the inter-stellar medium (ISM) by stellar 
evolution. 7/S(t*) is thus a monotonic function which grows 
from to 1 as t+ runs from the time galaxy formation starts 
to the present, and maps what fraction of was assembled 
up to a given lookback time. 

To track the evolution of the stellar metallicity we com- 
pute the total mass in metals inside stars as a function of t* 
and divide it by the total mass inside stars at the same time. 
This yields the time-dependent mean stellar metallicity: 

Z*(t*) = M(t*)Z*.i (5) 

t*,j>t* 

The computation of takes into account that the total 

mass inside stars changes with time, both due to the SFH 



1 Since our main goal is to compare the intrinsic evolution of 
galaxies, throughout this letter we consider ages in the context 
of stellar-evolution alone, ie., as if all galaxies were observed at 

2 = 0, such that t* is both a stellar-age and a lookback time. 



and to the time dependence of the returned mass. Accord- 
ingly, only populations older than t* enter the definition of 
Hi {t*), since younger ones had not formed yet, ie., fij(t+) = 
for t*j < t+. 



4 RESULTS 

The explicit time-dependence in equations Q and © pro- 
vides a description of the star formation and chemical histo- 
ries of galaxies which goes well beyond that obtained with 
the mere first moments of the age and metallicity distribu- 
tions. Given the numerous degeneracies which affect popu- 
lation synthesis of galaxy spectra one might be skeptical as 
to how useful they actually are. Simulations addressing this 
issue will be presented elsewhere. In what follows we report 
what is obtained in practice. 

4.1 Galaxy evolution in Z ne b-hins 

Our strategy is to chop the SF-wing into bins in Z ne b and 
study the evolution of galaxies from the statistics of rfi(t*) 
and Zi,{ti,~) in each bin. The fact that Z ne b correlates strongly 
with several physical and observational properties guaran- 
tees that galaxies inside each Z ne b bin are intrinsically sim- 
ilar, which in turn implies that statistics over such bins is 
meaningful. 

Fig |2K shows the mean mass-assembly functions ob- 
tained for each of the six Z ne b bins illustrated in Fig^s. The 
systematics is evident to the eye. Galaxies with low Z ne t, 
are slower in assembling their stars than those with high 
Z ne b- This behavior reflects the fact that, although most of 
the mass is assembled at early times (t*>l Gyr), galaxies 
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Figure 2. (a) Average mass assembly (r]+) and (b) chemical evo- 
lution (Z*) of stars in SF galaxies for the same six Z ne i,-b'ms 
defined in Fig Q). Each line represents a ti-by-t* average over 
all galaxies in the bin. A mild gaussian smoothing with FWHM 
= 0.2 dex in logi* was applied to the mean »7$(t*) and Z*(t*) 
curves, just enough to smooth out discreteness effects associated 
to the 25 ages in the base. 



at the tip of the SF-wing have larger star formation rates 
(SFR) per unit mass in the recent past than those at the 
bottom of the wing. This trend can be quantified in terms 
the ratio of the mean SFR in the last 60 Myr to that over 
the galaxy's history, ie, the so called Scalo's b parameter (eg, 
Brinchmann et al. 2004) . We find b to decrease by a factor 
of ~ 20 from the lower to the higher Z ne b bins. (SFRs are 
discussed in A07, where it is also shown that our synthesis 
based SFRs agree very well with those derived from the Ha 
luminosity.) Alternatively, one may quantify this behavior 
by computing the time at which the mean r?5(i*) function 
reaches 0.75, ie., the time at which the mass involved in star- 
formation reaches 3/4 of the total. We find that T75 increases 
from ~ 0.8 Gyr at Z neb ~ 0.31 (bin A in FigOJj), to 2 Gyr 
at Z na h ~ 0.77Zq (bin D), and 8 Gyr for the highest Z na h 
bin. Given that Z ne b correlates strongly with (Zaritsky 
et al. 2004; Tremonti et al. 2004; Lee et al. 2006), this behav- 
ior is ultimately a signature of galaxy downsizing (Cowie 
etal. 1996). 

The chemical evolution of stars in SF galaxies is shown 
in Figl^t). Again, the main message is clear: Galaxies with 
a metal poor ISM are low are slow in their stellar chemi- 
cal enrichment. In the lower Z n< , b bins, only reaches 
the current stellar metallicity level (~ 1/3Z©) as recently as 
t* ~ 100 Myr. Galaxies with larger Z ne b, on the other hand, 
have flatter and systematically higher curves. These 
massive systems have essentially completed their mass as- 
sembly and chemical evolution long ago, so that recent star- 
forming activity has a negligible impact upon Z+ft*). 

The fact that bins of increasing Z ne b have increasingly 



Figure 3. As Fig |2] but binning SF galaxies by their stellar 
mass, using five 1 dex wide bins, centered at (from bottom to 
top) log M^/Mq = 7.5, 8.5, 9.5, 10.5 and 11.5, which contain 
583, 4750, 29056, 46847 and 3550 galaxies, respectively. 

higher present Z* values confirms that nebular and stel- 
lar metallicities are related (SEAGal I; Gallazzi etal. 2005; 
A07). We nevertheless warn the reader not to overinterpret 
the quantitative relation between these two quantities, given 
that they are derived by completely different and hardly 
comparable means. In this paper we use [OIII]/[NII] as an 
indicator of Z ne b (equation^ , because of its evident relation 
to the position of the points in the BPT diagram. However, 
as emphasized by Stasihska (2006), strong line methods for 
deriving nebular abundances have been calibrated using gi- 
ant HII regions and their relevance for integrated spectra 
of galaxies has still to be examined, because of the role of 
metallicity gradients and of the diffuse ionized medium in 
galaxies. In addition, there is so far no study comparing neb- 
ular abundance calibrations with stellar abundances. There- 
fore, our Zneb values are not necessarily on the same scale as 
stellar abundances. Insofar as this paper is concerned, this 
is not a problem, since, irrespective of absolute values, it is 
the existence of a Z ne b-Zi, relation which is responsible for 
the ordering of the different chemical evolution curves in Fig 

4.2 Galaxy evolution in Mass-bins 

While Fig |5] illustrates the strong link between the stellar 
mass assembly and chemical evolution histories of galaxies 
and the present ionized ISM metal content, Z ne b is not the 
cause, but the product of galaxy evolution, so its use as 
an independent variable, although pedagogically useful, is 
questionable. If one had to chose a single quantity as the 
main driver of galaxy evolution this would surely be its mass 
(eg., Tinsley 1968; Larson 1974; Pagel 1997). 

We have thus repeated the analysis above binning galax- 
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ies by M* instead of Z ne b. Results are shown in Fig[3]for five 
f dex-wide mass bins centered at log M+/Mq = 7.5 . . . if. 5. 
Because of the flattening of the M* — Z ne b relation at high 
masses and the substantial scatter around the relation, par- 
ticularly at low M* and Z ne b, binning in M+ or Z n eb samples 
somewhat different populations of galaxies. On the whole, 
however, this alternative grouping scheme leads to the same 
general scenario outlined above, with low M* galaxies be- 
ing slower in their mass and chemical build up than more 
massive ones. A caveat is that M* is not a perfect tracer of 
the depth of the potential well, since it does not account for 
the mass in gas nor dark matter. In any case, the signature 
of galaxy downsizing is evident in Fig [3] massive galaxies 
today formed most of their stars and metals long time ago, 
whereas low-mass galaxies are still forming stars and metals 
actively. 

The signal of chemical evolution is present in all curves 
in Figs |3) and |3). Curiously, however, this evolution does 
not seem to start from very low Z* for massive, high Z„ e j 
galaxies (eg, bins E and F). We interpret this apparent lack 
of old, metal poor populations as due to a combination of 
lack of age resolution at high t* and an intrinsically fast evo- 
lution ocurring at early times. In low mass systems (eg, bin 
A), Z+(ti,) starts is At* ~ 2 Gyr long rise period at t* ~ 2-3 
Gyr, such that At*/£* is of order unity. The same process 
happening at t* > 10 Gyr would yield At*/t* < 0.2, and 
consequently be much harder to resolve. As shown above 
above, by t* = 10 Gyr massive galaxies had already assem- 
bled most of their mass, so their chemical evolution must 
also have occurred early on. This mixture of populations of 
widely different metallicities but similar ages results in Z± 
values at t* > 10 Gyr which are dominated by the metal rich 
stars already present at those times. Such an early and rapid 
formation period should lead to enhancement of elements 
synthesized in type II SNe (Worthey, Faber & Gonzalez 
1992; Thomas etal. 2005). Though we cannot map elemen- 
tal abundances with our method, we do detect a signal of 
increasing a enhancement towards the bottom of the BPT 
diagram by means of systematically increasing residuals in 
spectral bands due to a elements like Mg (A07), indirectly 
supporting the scenario outlined above. 



5 CONCLUDING REMARKS 

Our analysis shows that the present ISM abundance is in- 
timately connected to the past star-formation and chemical 
enrichment history of a galaxy. Whereas this is intuitively 
obvious from a physical point of view, it is not obvious at 
all that one should be able to recover this behavior from 
population synthesis analysis of integrated galaxy spectra 
of SDSS-like quality. We emphasize that, unlike previous 
studies such as that by Bica (1988), our fits impose no a 
priori chemical evolution constraints whatsoever upon the 
mixture of t*j 's and Z* j3 -'s in the base. Given the absence of 
constraints and the fact that the fits are completely indepen- 
dent of the emission line data used to compute Z ne b, the or- 
ganization of galaxies with different Z ne b onto the clear and 
systematically distinct stellar chemical enrichment patterns 
seen in Fig|3> is remarkable, revealing that fossil methods 
based on integrated spectra have reached a level of matu- 
rity which would be unthinkable just a few years ago (eg, 



Cid Fernandes etal. 2001). By taking the spectral synthesis 
output one step further than the trivial first moments of the 
i* and Z+ distributions, we have been able to uncover the 
chemical evolution of galaxies with an unprecedented level 
of detail for such a large and varied sample. This significant 
achievement was only possible due to the fabulous statistics 
of the SDSS, combined with state-of-the-art evolutionary 
synthesis models and reliable SFH recovery techniques. The 
sheer number of galaxies allows us to recover an unequivocal 
chemical evolution signal which would be at best doubtful 
for individual galaxies or even small samples. Combined with 
our estimates of the star-formation histories, the empirically 
derived chemical evolution patterns should provide valuable 
constraints for galaxy evolution models. 
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